CALT-68-2102 
hep-ph/9703213 
DOE RESEARCH AND 
DEVELOPMENT REPORT 



Model independent results for B — > Di(2420)£v and 
B -> D* 2 (2Am)iv at order AqcdM^ 



Adam K. Leibovich, Zoltan Ligeti, Iain W. Stewart, Mark B. Wise 
California Institute of Technology, Pasadena, CA 91125 



Abstract 



Exclusive semileptonic B decays into D\ and mesons are investigated 
including order Aqcd / f^c,b corrections using the heavy quark effective theory. 
At zero recoil, the Aqcd/^c,;; corrections can be written in terms of the 
leading Isgur-Wise function for these transitions, r, and known meson mass 
splittings. We obtain an almost model independent prediction for the shape 
of the spectrum near zero recoil, including order Aqcd/^c,& corrections. We 
determine r(l) from the measured B — > D\lv branching ratio. Implications 
for B decay sum rules are discussed. 
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The use of heavy quark symmetry |IJ resulted in a dramatic improvement in our under- 
standing of the spectroscopy and exclusive semileptonic decays of mesons containing a single 
heavy quark. In the infinite mass limit, the spin and parity of the heavy quark and that of the 
light degrees of freedom are separately conserved. Light degrees of freedom with quantum 
numbers sj l yield a doublet of meson states with total angular momentum J — si ± ~ and 
parity P = 7Tj. All semileptonic decay form factors of B mesons into either member of such 
a heavy quark spin symmetry doublet are given by just one function of w = v ■ v '. Here v is 
the four- velocity of the B and v' is that of the charmed meson. Moreover, for the B — > 
ground state to ground state transitions (these states have 1 = \ ), this universal function 
is normalized to unity at zero recoil 10-^]. Corrections to these model independent pre- 
dictions, suppressed by powers of AQ CD /m Cjfe , can be systematically investigated using the 
heavy quark effective theory (HQET) ||. 

In this letter we discuss semileptonic B meson decays into excited charmed mesons. Sur- 
prisingly, we find model independent predictions that hold even including order Aq C d / ^c,fe 
corrections. We concentrate on the doublet corresponding to s\ l = | + , which contains the 
£>i(2420) and the D^(2A60) mesons with widths around 20MeV. States in the sj l = \ + 
doublet can decay into D^n in an s-wave, and so they should be much broader than the 
D\ and D\ which can only decay in a <i-wave. (An s-wave decay amplitude for the D\ is for- 
bidden by heavy quark spin symmetry 0.) B — > D\iv and B — » D^lv account for sizable 
fractions of semileptonic B decays |0HH > an d are probably the only three-body semileptonic 
B decays, other than B — > £u, whose differential decay distributions will be precisely 
measured. 

The measured masses of various meson states containing a bottom or charm quark already 
give important information on HQET matrix elements. The D\ — Di mass splitting is only 
37MeV, suggesting that for sj l = | + states matrix elements involving the chromomagnetic 
operator are smaller than for the ground state (mp* — rriD = 140 Me V). The parameters A 
and Ai for the ground state multiplet can be related to the analogous parameters for D\ and 
D^i which we denote by A' and X[. We define the spin averaged masses % = (3m£>* +itid) /4 
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and m' D = (5rriD* + Stud^/S, and similarly for analogous mesons containing a bottom quark. 
Using the measured hadron masses |I0| , and identifying m' B = 5.70 GeV as the mass of the 
P}(5732) state, we find 

X[ — Ai = 2m c mb (m' B — TUb — fn' D + Id) / (mi, — m c ) ~ —0.34 GeV 2 , 
A' - A = m' D - fn D + (Ai - Ai)/(2m c ) ~ 0.35 GeV . (1) 

The matrix elements of the vector and axial currents {V^ = c^b and = 07^756) 
between B mesons and D\ or D\ mesons can be parametrized as 

(Di(v', e)| V" \B(v)) = y/m~^JnB~ [f Vl e*» + {f v ^ + f v ,v») (e* ■ v)} , 
(D 1 (v',e)\A^ \B(v)) = y/m Dl m B if A £^e* a v p v' y , 

(D* 2 (v', e) I A» \B{v)) = Jm Dl m B [k Al ^ a v a + (k A X + k Az v'») e* ap v a v?} , 

(D*(v', e)\ V» \B{v)) = ^m D *m B i k v ^e^>< • (2) 



Here and ki are functions of w. The differential decay rates for B — > D\ i v and B — > D\lv 
decays in terms of these form factors are, respectively (77 = uij) 1 j ttl b and r 2 = mD'/mg) 

dT a _ G 2 F \V cb \ 2 m%r\ 



dw 



+ 



48tt 2 



{w-r l )f Vl + {w 2 -l){f Vz +r 1 f v ,) 



dr 2 G 2 F \V cb \ 2 m%r 3 2 . 2 



dw 



144vr 3 



vV^l \2(l-2wr 1 +r 2 1 ) f 2 +(w 2 -l)l 



ki + (W 2 - 1) k 2 y 



[W 



l) 3/2 {3(l-2wr 2 + 



+2 



(w - r 2 ) k Al + (w 2 - 1) (k A3 + r 2 k A A 



(3a) 



(3b) 



The form factors and fcj can be parametrized by a set of Isgur-Wise functions at each 
order in Aqcb / m c,b- It is simplest to calculate the matrix elements in Eq. ([2]) using the trace 



formalism flll| . The fields P v and P* M that destroy members of the sj l = | doublet with 



four-velocity t> are in the 4x4 matrix 



l + i r 



(4) 



while for = | + the fields P" and P„* Miy are in 



(5) 



The matrices H and F satisfy jH v = H v = -H v f, j>F$ = F£ = -Ftf, F^ = 0, and 
r*F£ = 0. 

To leading order in AQ CD /m c b and a s 



cVb = h$T = r Tr [v a F% V H v ) 



(6) 



for matrix elements between the states destroyed by the fields in H v and F",. Here r is a 
dimensionless function of w, and is the heavy quark field in the effective theory. This 
matrix element vanishes at zero recoil for any Dirac structure T and for any value of r(l), 
since the B meson and the (Di, D%) mesons are in different heavy quark spin symmetry 
multiplets, and the current at zero recoil is related to the conserved charges of the spin- 
flavor symmetry. Eq. @ leads to the m Cj j — > oo predictions for the form factors /, and fc. 



given in Ref. [12 



At order AqcdA^c,& there are corrections originating from the matching of the b — > c 
flavor changing current onto the effective theory, and from order Aqcd/^c,6 corrections to 
the HQET Lagrangian. To leading order in a s , the current cTb is represented in HQET by 



,n,..«?fr-^?r + 5 i.r* + ...)e 



(7) 



For matrix elements between the states destroyed by the fields in F°, and H v , the new order 
Aqcd/^ c ,6 operators in Eq. (|7|) are 



h$ iDx T h® = Tr {SH F° T H v ) 
b$ V ifi x = Tr [S^l F$ T H v ) 



(8) 



Unlike B — > decays, since the .Dg) mesons and the B are in different multiplets, 
there is no relation between and <S^. The most general form for these quantities is 



S. 



(Q) 



A Q) v x + r^v' x + r^ lx \ + rP 9aX 



(9) 
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The functions r, depend on w, and have mass dimension one. They are not all independent. 
The equation of motion for the heavy quarks, (v ■ D) h)®' = 0, implies 



-(c) 



(<0 



WTi + r 2 v - r ; 







(10) 



Using T hff) = {Av v — A'v' u ) hffl T valid for matrix elements between the states 

in F£ and in H v , together with the equation of motion for the heavy quarks, and the 
constraints in Eq. (|10D, we obtain 



[w 



[W 



At zero recoil, rf = -r 4 w (l) = (A' - A) r(l) 



(b) 



.0>h 



T, 



(o) 



n 



(b) 



(wA' - A) t . 
(wA -A')t. 



(11) 



Next we consider the terms originating from the order AqcdA^c^ corrections to the 
HQET Lagrangian, 

1 



(12) 



2m Q v \>- ' ■ 2 

These corrections modify the heavy meson states compared to their infinite heavy quark 
mass limit. For example, they cause the mixing of the D\ with the J p = 1 + member of the 



i + 
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doublet. (This is a very small effect, since the D\ is not any broader than the D^.) 



The kinetic energy operator does not violate heavy quark spin symmetry, and therefore in 
B — > D\^D\ semileptonic decays its effect can be absorbed into a redefinition of r, which 
is used hereafter. For matrix elements between the states destroyed by the fields in F£, and 
H v , the time ordered products of the chromomagnetic term in 5C with the leading order 
currents are 



1+f 



(*) [f& r fcf )] (0) } = Tr n% F V Z io>* ±±L V H„ 



ijd*xT{ | a a , G af3 h®] (x) [T& T /#>] (0) } = Tr F°, T ±±L W#H V ) . (13) 

The most general parametrizations of IZ^^ are 



T>( c ) l c ) I ( c ) I ( c ) 

T>(b) (b) , (6) / , (6) / /-, A \ 

Only the part of antisymmetric in a and (3 contributes, when inserted into Eq. (p~3|) . 

The functions rji depend on w, and have mass dimension one. (Note that g aa 1i3 is dependent 
on the tensor structures included in Eq. flU]).) All contributions arising from the time ordered 
products in Eq. ([13]) vanish at zero recoil, since v a F£ = 0, and v a {\ + ^)a a/3 (l + f) = Of] 

Using Eqs. (ffl)-(O), it is straightforward to express the form factors fi and ki parametriz- 
ing B — > D\lv and B — > D\lv semileptonic decays in terms of Isgur-Wise func- 
tions. We use the constraints in Eqs. ([UJ) and (|TTf) to eliminate r 3 and r 4 . The 
form factors in Eq. (Q) depend on and rff 1 only through the linear combinations 
n = (2w + 1) rf } + r 2 (6) + 2(A' - wA) r and r} b = 6r}[ b) - 2(w - 1) rf^ + 77^. Denoting 
e Cj b = l/(2m cfc ) and dropping the superscript on and 77^, we obtain for the 5 — ► Lq Iv 
form factors 

v^/a = -(«; + l)r - e 6 [(w - l)r 6 + (w + 1)^] 

+ £c [3(w; - l)(n - r 2 ) + (w + 1)(2^ + 3r/ 3 ) - 4(wA' - A)r], 

V6f Vl = (1 - w 2 )r - £ b (w 2 - l)(r 6 + Vb ) 

+e c [{w 2 - 1)(3t! - 3r 2 + 2^ + 3%) - 4(iu + l)(wA' - A)r], (15) 

V^/va = -3r - 3£ b (r b + - s c [(Aw - 1)t x + 5r 2 + IO771 + A(w - 1)772 - 5r/ 3 ], 

\fcfv 3 = (w - 2)r + e b [(2 + w)r b - (2 - w)i] b } + e c [(2 + iy)n + (2 + 3u>)r 2 
-2(6 + w)rn - A(w - 1)773 - (3iu - 2)773 + 4(wA' - A)r], 

The analogous formulae for 5 — > D* 2 iv are 



* Order Aqcd/^c corrections were also analyzed in Ref. [13|. We find that T4 (denoted £4 in |lj 
does contribute in Eq. @ for T = 7Ar, and corrections to the Lagrangian are parametrized by 
more functions than in 
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ky — t e b (r b + r] b ) - e c (ri - t 2 - 2?7i + 773), 
k Al = -(1 + w)t - e b [(w - l)r b + (1 + w)i] b ] 

-e c [(w - l)(n - r 2 ) - (w + l)(2?7i - rj 3 )}, 
k A2 = -2e c (r 1 +r ]2 ), (16) 
k Aa = t + e b (r b + r] b ) - £ c (t! + r 2 + 2?7i - 2?7 2 - 773). 

The allowed kinematic range for B ^* D\iv decay is 1 < w < 1.32, while for B — > D\iv 
decay it is 1 < w < 1.31. Since these ranges are fairly small, it is useful to expand the 
differential decay rates in Eq. (Q) simultaneously in powers of w — 1 and Aqcd/^c,6- 

12 IT/ 12 ^,5 ^,3 



dl\ 2 G 2 F \V cb \ 2 m%rl 2 _ 

— — ~ — — Vw 2 - 1 

dw 487T 3 



(17) 



x 



In a; 1*2 (2 — 0, 1, 2), we keep terms up to order (AQ CD /m Cjfe ) 2_ \ Eqs. (|), (fl5|), and ([IB]) yield 

4 0) = 32£ c 2 (A'-A) 2 (l-r 1 ) 2 , 
X W = (8/3) [(1 - n) a + 4e c (A' - A) (3 - 4 n + r 2 ) 
- 2(1 - ri) 2 {2e c rn + 3e c r] 3 - e b 7] b )/r] , 
xf ) = (8/3) [(3 - 4n + 3r 2 ) + 2(1 - n) 2 r'/r] , 

4 0) = 0, (18) 

= (40/3) (1 - r 2 ) 2 [1 - 2(2£ c 7?i - e cV3 - e 6 %)/r] , 
x { ? = 8 (3 - 8r 2 + 3r 2 ) + (80/3) (1 - r 2 ) 2 r'/r . 



Here r, r' = dr/dw, and 77, are evaluated at w — 1. The values of r/j that occur in Eq. (|T8"|) 
are not known. Since the D\ — Di mass splitting is very small, and model calculations 
indicate that the analogous functions parametrizing time ordered products of the chromo- 
magnetic operator for B — > iv decays are tiny [[14]], hereafter we neglect the corrections 
parametrized by rji. 

The value of r(l) can be determined from the experimental measurement of the 



B — > Di £ v branching ratio. We use the average of the ALEPH and CLEO || results, 
Br(fi -+D x lv) = (6.1 ± 1.1) x 1(T 3 , to obtain 

r(l) = 0.55 ±0.05. (19) 

(We used A' - A = 0.35 GeV from Eq. flip, r B = 1.6 ps, \V cb \ = 0.04, and m c = 1.4 GeV.) 
To get Eq. (||) we assumed r'(l)/V(l) = —0.8. r(l) has little sensitivity to this choice. 
Allowing -1.2 < t'(1)/t(1) < -0.5 only affects the central value in Eq. (H) by ±0.01. The 



ISGW nonrelativistic constituent quark model predicts r(l) = 0.54 in surprising agreement 
with Eq. (H) Hp]. 



The ALEPH and CLEO analyses assume that B — > Di i v X is dominated by B — > Di i u, 
and that Di decays only into D*ir. If the first assumption turns out to be false then r(l) 
will decrease, if the second assumption is false then r(l) will increase compared to Eq. (|i~9|). 

Even though £ C (A' — A) ~ 0.12 is small, the term proportional to it comes with a large 
coefficient, and dominates the value of x± . Numerically, this Aqcd/^c correction to x^ 3 
is 1.8, while the part that survives in the m Ci j — > oo limit is 0.8. Note that the part of the 
Aq CD /m^ correction to Xi that involves A', A, and r'(l) is only 0.27±0.10 for the previously 
mentioned range of r'(l) (using A = 0.4 GeV |T(|). The order Aqcd/^c terms that involve A' 



and A change x\ by less than a third of its leading order value for —1.2 < r (l)/r(l) < —0.5. 

After Eq. (|I"2|) , we absorbed into r the form factor that parametrizes time ordered prod- 
ucts of the kinetic energy operator with the leading order currents. While A' x is quite large 
(see Eq. (|])), this is probably a consequence of the Di and D\ being p- waves in the quark 
model, and does not necessarily imply that the heavy quark kinetic energies significantly 
distort the overlap of wave functions that yield the form factors. If we had explicitly in- 
cluded the time ordered product involving the charm quark kinetic energy, the leading term 
in Xi would change from (1 — r^) 2 to (1 — r x ) 2 (1 ± 2r]^ c e^). Even taking r]^ c = ±A' changes 
the extracted value of r(l) by less than 0.04. So this Aqcd/^c correction to is likely to 
be much smaller than the term proportional to e c (A' — A) explicitly shown in Eq. flIB|). It 
is important to have experimental data on the w-spectrum of B — > D\ i v decay to test the 
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hypothesis that the effects of the kinetic energy and AQ CD /m 2 corrections are not large. 

The order Aqcd/^c^ corrections calculated in this letter are also important for the 
prediction of R = Bi(B -> D%£p)/Br(B -> D x lv). As R is sensitive to r / (l)/r(l), we 
shall explicitly display the dependence on r'(l). In the m c ^ — > oo limit, expanding to 
linear order in r'(l)/r(l), we obtain R = 1.89 + 0.51 t'(1)/t(1). Including the AQ CD /m 2 



correction to , and the Aqcd/^c correction to , and expanding again to linear order 
in r'(l)/r(l), yields R = 0.79 + 0.30 r'(l)/r(l). This suppression of R compared to the 
infinite mass limit is consistent with experimental data. (It is possible that part of the 
reason for Bi(B -^D*£i?X)x Br(L>* -> D^ir) < (1.5 - 2.0) x 10~ 3 is a suppression of 
Bi(D* -> DWtt) compared to Br(L>i -> D*tc).) 

Our results are important for sum rules that relate inclusive -B — > X c I v decays to the sum 
of exclusive channels. The Bjorken sum rule [17ll2| for the slope of the B — > iv Isgur- 



Wise function becomes p 2 = —d£/dw\ w= i > 0.25 + r 2 (l) = 0.55. Note that 2r 2 (l)/3 arises 
from the D 1 , D* 2 doublet, while r 2 (l)/3 is due to the broad sj l = \ + doublet, (Dq, D{). 
(Using the equality of the leading Isgur-Wise functions for these multiplets in the quark 
model, valid for any spin-orbit independent potential.) A class of zero recoil sum rules were 



considered in Ref. ||18||. The axial sum rule, which bounds the B — > D* form factor that 



measures \V c b\, receives no corrections from either the | + or the | + doublets. The J p = 1 + 
states contribute to the vector sum rule, which bounds the Ai parameter of HQET. This 
bound is strongest in the limit m c ^> m b ^> Aq C d, where the D\ state does not contribute. 

Q _|_ -| | 

The equality of Isgur-Wise functions for the | and | doublets in the quark model implies 
that B -> D\ transition modifies the bound to \ x < -3A 2 - 3(A* - A) 2 r 2 (l) ~ -3A 2 - 0.09. 



(Here A* ~ A+0.31 GeV [19j is the analogue of A for the \ + states.) Perturbative corrections 



to these bounds can be found in |20|. 



In this letter we analyzed B — ► Di(2A20) £ i? and B — > 0^(2460) iv decay form factors 
at order Aqcd/^c^- At zero recoil, all Aqcd/^c,6 corrections can be written in terms of 
the m Cj b — > oo Isgur-Wise function for these transitions, and known meson mass splittings. 
With some model dependent assumptions, we predicted the shape of the spectrum near 



zero recoil, including order Aqcd/^c,?) corrections. Testing these predictions will constitute 
an interesting check on our understanding of exclusive semileptonic decays based on the 
HQET. Similar results hold for semileptonic B decay into the broad sj l = | + charmed 
meson multiplet, and for the semileptonic A& decays into excited charmed baryons. These 
will be presented in a separate publication. Perturbative QCD corrections and nonleptonic 
decays (using factorization) will also be considered there. 
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